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ABSTRACT 

, Context. The dynamics contained in magnetized layers of exoplanet atmospheres are important to understand in order 

to characterize what observational signatures they may provide for future observations. It is important to develop a 
' framework to begin studying and learning the physical processes possible under those conditions and what, if any, 

' features contained in them may be observed in future observation missions. 

Aims. The aims of this study is to formally derive, from scaling arguments, a manageable reduced set of equations for 
analysis, i.e. a magnetic formulation of the equations of quasigeostrophy appropriate for a multi-layer atmosphere. We 
check these derived equations for consistency with respect to similar equations currently used in the literature like the 
magnetized shallow water equations and their precursors. The main goal is to provide a simpler theoretical platform to 
explore the dynamics possible within confined magnetized layers of exoplanet atmospheres. 

Methods. We primarily use scaling arguments to derive the reduced equations of "magnetoquasigeostrophy" which 
assumes dynamics to take place in an atmospheric layer which is vertically thin compared to its horizontal scales. 
^ Furthermore, the derivation exploits the fact that the Rossby Numbers of the emergent flows are small and that the 

, Cowling Number is also near 1, the latter of which measures the relative content the energy per unit volume contained 

(— I ■ in the magnetic field to the corresponding kinetic energy of the flow. 

■ Results. The magnetized incarnation of the quasigeostrophic equations for a two-layer system are fully rederived from 

I ' scaling arguments. The resulting equation set retains features existing in standard shallow-water magnetohydrodynamic 

Q , equations but are absent in more classical derivations of the quasi-geostrophic limit, namely, the non-divergence of the 

?H , in-plane components of the magnetic field. We liken this non-divergence of the in-plane magnetic fields as indicative 

C/j ' of a quantity whose behaviour mimics a two-dimensional "pseudo" -magnetic monopole source. We also find, using the 

' same scaling argument procedures, appropriate limits of the fundamental parameters of the system which yield reduced 

equations describing the flow dynamics primarily characterized by magnetostrophic balance. 

Conclusions. The standard scaling arguments employed here show how traditional magnetized quasigeostrophic equations 
. connect to their magnetized shallow water forms. The equations derived are amenable to analysis using well-known 

' techniques. 

00 ' Key words. Hydrodynamics, Magnetohydrodynamics, Exoplanets 
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1. Introduction 

It is now generally accepted that hot extrasolar giant planets ( "hot- Jupiters" ) are commonplace in the Galaxy. Many 
orbit their parent stars in so close that recent calculations by Koskinen et al. (2010) indicate that the upper atmospheric 
layers of these hot- Jupiters are sufficiently irradiated by the parent star's UV radiation field (and therefore sufficiently 
ionized) to treat the gas as a plasma. It is important therefore to understand what sorts of observational indications 
these upper layers may show to future telescope missions. Model atmosphere calculations show that these upper layers 
I are likely to be stably stratified (like the Earth's stratosphere) and, therefore, it is justifiable to use global circulation 
d • models to describe and explore their flow dynamics (Cooper & Showman 2002, Cho et al. 2003, Menou & Rauscher 
2009, to name only a few). Preliminary examinations of the response of magnetized shallow- water models of exoplanets 
reported by Cho (2008) indicate that under not-too unreasonable conditions, the atmospheres of irradiated exoplanets 
could settle into a collection of stationary, planetary scale, magnetized vortices with zero total circulation. If this turns 
out to be a robust feature for hot irradiated Jupiters (for instance) then there will be obsevational consequences which 
may be detectable by the next generation of space missions. 

The aim of this work is to begin to formulate a systematic approach toward understanding the dynamics of stably 
stratified flows on exoplanets under the influence of magnetohydrodynamic effects which are usually absent in standard 
meteorological modeling. The simplest place to begin exploring these kind of dynamics is in a quasi-geostrophic ("QG" 
for short) framework which is a setting representing dynamics (i) occurring on synoptic scales (horizontal scale dynamics 
which are much smaller than the planetary radius), (ii) described by order 1 Burger numbers and (iii) characterized by 
small Rossby number flows. The last of these is a measure of the ratio of the typical rotation time of the planet to the 
circulation time of a synoptic scale vortical structure. The Burger number is also a ratio between the planet's rotation 
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time to the gravity wave propagation time across the synoptic scale. If such an atmosphere is also magnetic, then another 
important number governing the quality of the dynamics will be the Cowling number which measures the the relative 
importance of the Lorentz force. The Cowling Number here will be understood as the ratio of the energy contained in 
the magnetic field to the kinetic energy in vortical motions. 

When the Cowling number becomes an order 1 quantity QG dynamics will be substantially modified. Examination 
of magnetohydrodynamic modified quasigeostrophic flow goes back to the original series of studies in the thesis work 
of Gilman (Gilman 1967a-c) motivated by the problem of solar differential rotation. The model equations derived in 
Gilman's original studies are the magnetohydrodynamic analog of the classical quasigeostrophic equations on a /3-plane. 
In these equations the magnetic field is strictly two-dimensional and horizontal, where it is assumed that the magnetic 
energy content contained in this layer is outweighed by the energy contained in the hydrostatic configuration. As a 
result, the in-plane magnetic field does not figure into the lowest order geostrophic balance. Furthermore, the field lines, 
although they behave two-dimensionally, can generate vertical field lines through the action of vertical layer undulations. 
The variation scales of these vertical field fiuctuations are, however, small compared to the horizontal divergences (at 
next order) and, thus, though vertical field lines may be generated they play no dynamical role in the equations of motion 
at lowest order. 

These equations also support Rossby waves and Alfven waves, often referred to by their hybridized forms as 
'hydromagnetic-planetary' waves (Acheson & Hide 1973). Unlike the situation in classical QG, however, the magneto- 
quasigeostrophic equation (MGQ) set cannot be characterized by a conserved quantity like the potential vorticity in 
QG. However, the magnetic potential function (sometimes called a "flux" -function) is a materially conserved quantity of 
these equations. Also, motivated by meteorological modeling, Gilman (1967b) extended the basic MGQ formalism onto 
a multi-layer "stacked" model where the density in each layer are constant possibly differing from each other. As it is 
one of our goals to formulate something similar, Gilman's study shall serve as a guiding procedure in this respect. 

The magnetohydrodynamic version of classical shallow water equations (hereafter SWMHD) was developed in Gilman 
(2000) in order to address the problem of the solar tachocline (Spiegel & Zahn, 1992). The SWMHD are model equations 
confined to a thin-shell of magnetized gas which is stably (or marginally) stratified. In the familiar non-magnetized 
setting, the shallow water equations retain the action of gravity waves while the more simplified QG models do not. 
Gravity waves are absent in the latter because the Burger numbers are order 1 in QG, and since the timescales for 
dynamics in QG are much longer than the local planetary rotation time (because the Rossby number is small), it means 
that gravity waves are filtered out. In this sense, one may regard the shallow water equations as containing "more" 
physics over the QG set. The same is true of the SWMHD model in that it retains the hydromagnetic planetary waves 
found in the MGQ set and it further supports magnetically modified gravity waves whose properties and structure has 
received recent scrutiny (Schecter et al. 2001, Zaqarashvili et al. 2007, Heng & Spitkovsky 2009). 

Besides the absence of magnetically modified gravity waves, the SWMHD equations and the original MGQ models also 
diverge from one another in content because the horizontal divergence of the magnetic fields are not zero in SWMHD as 
they are in the original classical MGQ equations. This means that, for example, the dynamical evolution of the horizontal 
fields cannot be formulated in terms of the evolution of a scalar 'fiux- function' as is usually done for horizontal-divergence 
free fields. It is shown in this study that one can view this non-divergence of the horizontal fields as corresponding to 
the presence of a globally conserved, two-dimensionally distributed "pseudo" -magnetic monopole charge, q^"^\x,y,t).^ 
Therefore, in this sense, one may characterize the solutions contained within the various model approximations according 
to whether or not (7'™^ is zero. For example, the MGQ developed in the original papers of Gilman from 1967-8 have 
= while the SWMHD models have q^™) 7^ in general. 

What is achieved in this work: Since classical QG offers a sound conceptual platform to focus on and understand the 
dynamics of (mainly) small Rossby number vortical flows, it would be logical to start studying magnetized exoplanet 
atmospheres from a similar vantage point. It makes sense to begin this exploration from a MGQ framework and, once 
having uncovered some of the dynamics that this set of equations contain, verifying that some facet of these dynamics are 
contained in the SWMHD equations. However, some ambiguity persists as to the fate of the so-called "pseudo" -magnetic 
monopole distribution in the reduced MGQ system as it is identically zero in Gilman's original studies but clearly present 
in SWMHD. Does the MGQ framework contain some features of this "pseudo" -magnetic monopole field g^™^? The answer 
is yes and it is one of our goals here to develop an intermediate framework to study magneto-vortical dynamics free of 
gravity waves but not necessarily free of the g*-™-* field. In this work the MGQ equations for arbitrary number of vertical 
layers of constant density is re-derived using familiar scale analysis techniques as, for instance, found in Vallis (2006) (see 
also Pedlosky 1987). We display the derived equations for a two-layer model. We show how the various aforementioned 
features between equation sets relate to one another. 

The procedure implemented to derive the MGQ equations are sufficiently general also to establish the scalings which 
lead to magnetostrophic balance (Acheson & Hide, 1973). Magnetostrophy is a limiting form for small Rossby numbers 
in which dominant dynamical balance exists between Coriolis, pressure and Lorentz forces. The magnetostrophic ap- 
proximation has been extensively explored to study the Earth's geodynamo and associated turbulence (e.g. Fearn 1997, 
Moffat 2008). In our analysis and overall setting we find that the magnetostrophic balance at low Rossby numbers is 
achieved when (i) the Cowling number is sufficiently large in proportion to the inverse Rossby number and (ii) the Burger 

^ This unfortunate tongue-in-cheek description is used for the lack of a better expression to describe how the non-zeroness 
of the horizontal divergence of the magnetic field would "appear" as a monopole charge source, albeit, in two-dimensions. Of 
course no such thing exists as the non-vanishing feature of the horizontal divergences of the magnetic fields means only that there 
correspondigly exists a vertical field with non-zero vertical gradient. 
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number is small in proportion to the Rossby number. After establishing this fundamental magnetostrophic balance we 
proceed to derive a set of self-contained equations describing dynamics which are analogous to the MGQ set. We only 
derive the resulting equations here and shall return to a more thorough discussion of them in a follow-up study. 

In Section 2 the equations and assumptions are presented. In Section 3 the equations of motion are expanded around 
a given latitude, scaled and analyzed. In this section we formally define the non-dimensional parameters of the system 
and scale the equations of motion. Most importantly, we highlight the fundamental relationships that must be met 
between the various parameters in order to attain the aforementioned MGQ and magnetostrophic balancing. The actual 
expansion procedure completing the derivation of the MGQ equations is detailed in Section 4 while the same is done for 
the magnetostrophic set in Section 5. Because the main focus of this study is about the MGQ equations, in Section 6 
they are summarized followed by a brief discussion demonstrating how a MGQ model reduction does exist which retains 
the g^™^ field and how, most importantly, the classical MGQ system is recovered when gf^™^ is set to zero. 



2. Equations and assumptions: an overview 

The general equations of motion of an ideal, incompressible MHD fluid in a frame rotating with rotation vector ft are 

4^ + 2f2curlU = igradP - gr + ijcurlB, (1) 
dT p ' p 

div-U = 0, (2) 

— + curl(U curlB) = 0, (3) 

divB = 0, (4) 

where "div" , "grad" and "curl" are the corresponding three dimensional operations of divergence, gradient and curl 
respectively. These equations will be considered in the next section in a Cartesionized representation around a point 
at some latitude located suitably away from the planet's equator. In the Cartesionized frame, X represents the zonal 
(longitudinal/east-west) direction while Y represents the meridional (latitudinal/north-south) direction while Z rep- 
resents the vertical (i.e. planetary radius) direction. The corresponding velocities are (respectively) U, V, W while the 
magnetic fields are Bx, By, Bz- P represents the pressure. We consider the dynamics of these equations subject to the 
following setting and restrictions: 

1. The planetary /3-plane, nominally centreed on a latitudinal zone which is significantly away from the planet's equatorial 
zone. 

2. The dynamics of small Rossby number disturbances (i.e. Ro <C 1) are of interest so that we may develop the equations 
as a power series in this quantity . 

3. As is standard practice in atmosphere modeling, we break the atmosphere up into several layers. We assume that the 
background density fields in each vertical layer is constant so that each layer is effectively incompressible as implied by 

The layers, however, may have differing densities each given by pi where i denotes the layer under consideration. 

4. We shall allow for either none, some or all of the layers to be electrically conducting so that MHD will be a good 
description of the dynamics in those layers that are magnetically active. This is why the induction equations for ideal 
MHD are included. 

In the majority of the analysis performed in this study the dynamics are treated in two-layers. However, for the 
sake of generality of the procedures we implement, let us consider the possibility of an arbitrary number of constant- 
density layers with, as yet, unassigned conductivity. Let us suppose that the heights delineating the transition from one 
constant density layer to another is given by T-Li < T-L2 respectively where Z = represents the nominal "bottom" of the 
atmosphere (see Figure 1). Then the prescription for the current J will be represented in this way by 

ipN, Hn-i < Z < Hn 

n,-i <z<u, (5) 

The functions Lpi will take on either the values or 1 depending upon the layer i. 

Sections 3-4 detail the derivation of the MGQ equations with explicit presentation of a two-layer model. The derivation 
follows the standard procedure outlined in Vallis (2006, § 5.3, pgs. 207-215) while exploiting some of the scaling argu- 
ments invoked by Gilman (1967a). We introduce here some of the relevant scalings and corresponding non-dimensional 
parameters appearing in the system analyzed. The planetary rotation as viewed from the latitude in question is scaled 
by If the typical meridional/zonal velocity scales are given by lA then we may define the Rossby number 



J = (^curlB; Lp ■ 



If the planetary rotation is given by SI and if one is at latitude A then the rotation normal to that latitude is Slo = O, sin A 
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where /o = 2ilo and R is the planetary radius and C represents the horizontal "synoptic" scale of the dynamics. 
Dynamic flow timescales are assumed to be 1/Ro times longer than the Coriolis timescale 1/fo - this is the basis of 
the quasigeostrophic analysis. The magnetic field strength is scaled by B. This then leads to a natural non-dimensional 
quantity 

1 

— 

also known as the Cowling number. The scale density p is equated to the density of the lower atmosphere layer. C is 
related to the inverse square of the usual j3 parameter frequently referred to in plasma physics studies. The Cowling 
number can be understood as a measure of the energy contained in the magnetic field versus that contained in kinetic 
motions. It can also be understood as the square of the magnetic Mach number, i.e. C = Ua^ /LP, since U\ = B'^/A-Kp is 
the Alfven speed squared. 

The atmosphere also has the following vertical length scales that are important: the overall vertical extent of the 
atmosphere %, and the vertical length scale of atmospheric /luctuations h. Therefore, the zonal and meridional dynamical 
lengths {X, Y) are scaled by C while the vertical scales [Z) are characterized by %. Given that the flow speeds are O (JA) 
it follows that their typical characteristic timescales (T) are Ro"-*^ times longer than I//0, consistent with the assertion 
made above. These are all formally reintroduced in the discussion appearing in the next section. Finally, the Burger 
number, Bu, measuring the relative importance of gravity, is given by the relationship 

Bu^- 



3. Scalings and analysis 

The local Cartesionization of the equations of motion around a latitude Aq is given by, 

^ + 20zxU = -Vn - ffz + -^(B • V)B, (6) 

dT p Attp 

V • U = 0, (7) 

^ = (B.V)U, (8) 

V-B = 0, (9) 

where the total pressure 11 is given by 

n = p + ^B2. (10) 

The coordinate Z coincides with the effective gravity direction. The component of gravity g will be taken to be constant. 
In this construction X, Y denote the zonal and meridional coordinates respectively. 

Before analyzing these equations on the Cartesion plane, we explicitly nondimensionalize all quantities appearing in 
the equations of motion according to the dimensional quantities stated above. Henceforth, all lower case Latin symbols 
(some with hats over them) denote the corresponding non-dimensional quantity. The primary parameter that will be 
used for the following expansions will be the Rossby number, Ro. In terrestrial and planetary studies the Rossby number 
is typically substantially less than one, i.e. Ro <^ 1. U is the typical velocity scale of the horizontal planetary scale 
atmospheric motions, that is to say, 

U^Uu, V ^ Uv, (11) 

where m, v are the non-dimensionalized order one representations of the longitudinal and latitudinal velocities. Since U 
measures an effective overturning speed of the synoptic scale vortex structures in the atmosphere, the corresponding time 
scale associated with it scales as Ro/f2o- Thus we scale the horizontal lengths and time according to: 

X ^jCx; Y ^ Cy, T ^ t/(RoOo) 

where x,y,t are now understood to be non-dimensionalized quantities. The vertical scale is measured by H and this will 
be considered small compared to C. In problems in which the atmosphere is broken up into subplaycrs "H represents the 
full vertical extent of the whole atmosphere under consideration. In order to recover the quasi-geostrophic ordering the 
smallness of the ratio H/jC needs only be ^ 1. Correspondingly, the vertical velocity will also be assumed to be small in 
proportion to the ratio H/C Thus we scale these two quantities accordingly as 

Z^Hz, W-^ Uw 
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where z and w are order 1 non-dimensional quantities. The variation of the Coriohs parameter (which is the usual function 
of latitude on a planet) is written as 



R 

where f2o is the evaluation of the projected planetary rotation at latitude Aq in which yC/R = A — Aq. ilo — f^sinAo 
and /3 = cos Aq, where Vl is the planetary rotation rate. The ratio £/i? will be formally shown to be O (Ro) in Section 4 
in order to retain the planetary-/? effect. The magnetic field strength will be measured by B. Thus, in keeping with the 
formalism developed by Gilman (2000), the horizontal fields are scaled by this, 

Bb^, By -> Bby, 

while, because of the small aspect ratios under consideration, we shall only consider vertical fields which are smaller than 
the horizontal components in proportion to the ratio H / 

H 

Bz —bz. 

The vertical extent of the atmosphere is subdivided into sublevels with vertical scale Hi where i labels the level number. 
For our considerations we assume that all AHi = Tii — Tii-i, representing the thicknesses of each layer, are an order 1 
fraction of the overall vertical scale H. We let hi represent perturbations about the mean level height Hio, i.e., 

Tii — Hio + /ij. 

If we assume that the deviations scale by a characteristic length h, then we may introduce another parameter 5 which 
measures the the deviations against the vertical scale of the atmosphere 

^4- (12) 

In typical quasigeostrophic scalings this factor 5 is assumed small, usually O (Ro) for Ro small. 



3.1. Analysis of the vertical momentum equation 

We seek to develop reductions of the equations of motion which are dynamically in hydrostatic equilibrium. This means 
assuming the pressure scaling to be 

n ^ .gT^pII, 

where p is the density scale, typically of the lowest atmosphere layer. Examining the vertical component of the momentum 
equation with the scalings thus far presented reveals 

(^* + "-^ + "^^-)^ = Bu'^ (^a.n + lj +Ro2c^ (^^^ {h-v + bzdz)bz. (13) 

Where the Burger number, Bu, and the Cowling number, C, are as defined at the end of Section 2. To insure that 
hydrostatic balance is always dominant we shall assume Bu always dominates the other terms appearing in (I13p . This 
means to say specifically that unless 

Bu^ > Ro^ , and Bu' > Ro' C, (14) 

the assumption of hydrostatic balance will fail. For planetary atmospheres this criterion is generally met since H/ C <^ 1. 
The new piece here is the second disparity involving the Cowling number. We return to this shortly below. As of this 
point we have said nothing about the specific orderings required of the both the Burger and Cowling numbers. In general 
for this study we shall assume that Ro and the ratio H/C are sufficiently smaller than 1. It follows that the lowest order 
balance of the vertical momentum equation (jl3p is hydrostatic, 

5^11+^=0. (15) 
P 

Henceforth we shall treat p as being non-dimensionalized by p which means, in practice, that p/p will be replaced 
everywhere by p. When multiple layers are handled we shall characterize the density as being constant in each layer 
with value pi. Furthermore, to be consistent with what was stated earlier, the density of the lowest layer will always be 
unity, thus pi = 1. The procedure performed here essentially mimics the strategy employed in Vallis (2006). We allow 
the top layer to move about freely and assume the pressure is zero on the top surface. The solution for the pressure 
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for an N-layered atmosphere (with a flat bottom located at z = 0) expressed in terms of the non-dimcnsionalized lid 
coordinates Hi is given by 



J2k=N Pki^k — Hk-l) + PN-i{Hn-1 — z), Hn-I > Z > 



Pn{Hn - z), 



Hn > z>H, 



N-l 



n / -dz = { 



Ef=j+i Pk{Hk - Hk-i) + Pi{Hi - z), 



H, > Z> H; 



(16) 



Y^k=2 Pk{Hk - Hk-i) + pi{Hi - z), 



iJi > z > 0. 



Note that in the above formulation the summation procedure begins with layer i = N — 1 and that it is explicitly absent 
in the top, i = N, layer. As was motivated earlier, the levels Hi are further written in the form 



H, = H 



iO 



Shi, 



(17) 



where the set {Hio} are constants. Henceforth the procedure outlined in the remainder of this study will focus only on 
a two-layer system with the understanding that it straightforwardly generalizes to a multi-layer configuration. Given 
the functional form postulated above we find that horizontal gradients (V = Tidx + ydy) of the pressure are given by 



Vn = 5Vn = (5V< 



H2> z>Hi 



III = + (pi - P2)^i, Hi> z>0. 



(18) 



3.2. Analysis of the horizontal momentum equations 

Now we are prepared to move onto the horizontal momentum equations. With the scalings proposed above the equations 
become 

Ro^ {dt + u-V + wd^) u + 2Ro (^1 + ^Py^ z x u = 

-(JBu^ivn + Ro^-(h ■ V + b,d,)h + o( Ro^- ) , (19) 
p p \ RJ 

in which u = uic + vy and b = bx'x. + byy. We note that the curvature terms neglected at this stage come in at an ordering 
of O(Ro^) -OiC/R), and since £/R will be assumed to be <C 1, the curvature terms will be neglected throughout. It is 
at this stage that we explore various balances: 



3.2.1. Quasigeostrophic/magnetoquasigeostrophic scalings 

We assume that the deviations from the mean height scale to the Rossby number, i.e. 6 = O (Ro) . Although it would 
be enough to require that C/R ^ 1, in order to include the planetary /3 we require that C/R O (Ro) . In Section H] 
we detail the full procedure leading to the quasigeostrophic reduction, however, we show here the leading order balances 
to illustrate the main flavour of the results. Expanding the pressure and velocities within each layer i in powers of Ro 

according to, fl^ = flf^ + Roflf ^ + ■■■, = uf^ + Rouf ^ + ■ • • and b, = bf ' + Rob|^^ + • • •, then it follows that the 
leading order terms of which are ©(Ro), reduce to the fundamental statement of geostrophic balance, that is to 
say, 

2z X uf ^ = -Bu^lvnf \ (20) 

Pi 

Because we consider the density in each layer to be constant, taking the curl of the above equation immediately reveals 
that the horizontal divergence of the lowest order velocity field is identically zero 

) + = 0. 

Therefore it follows from this that the lowest order vertical velocity wf^'' must zero. Requiring S to be less than order 1 
actually follows from the horizontal divergence free condition of the lowest order horizontal velocity field. This consistency 
is formally shown in Section |4l 

In order to understand how magnetic effects enter the resulting equations we carry out the expansion to order Ro^ 
to find, 

{d, + uf) . V)uf ' + 2z X + I3yz x ' = -Bu^lvftf ^ -I- ^(b|") ■ V)bf 

Pi Pi 

^ To avoid any ambiguity: the general form of the disturbance pressure would be given by the general formulae tlf''' — pihi + 
J2k^i+i Pk{hk — hk-i) with n^' = pnhn. 
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Table 1. Summary of scalings and lowest order balances for hydrostatic disturbances. All balances shown assume Ro 
■C 1. The quoted extreme values in describing magnetostrophic balance is the same as saying the Acheson and modified 
Burger numbers (denoted in the text by A and Bu respectively) are 0(1). 



Balance Name 


C»(Bu) 


0{6) 


o{c) 




Lowest Order Balance 


QG 


1 


Ro 


< 1 


Ro 


Coriolis, Pressure Gradient 


MQG 


1 


Ro 


1 


Ro 


Coriolis, Pressure Gradient 


Magnetostrophy 


Roi/2 


1 


Ro"^ 


< 1 


Coriolis, Pressure Gradient, Lorentz Force 



The above is a restatement of Eq. p7p where in arriving at this equation it is argued that the horizontal components 
of the magnetic field are independent of z at lowest order. All but the last term on the RHS of the above expression 
constitute the next order corrections which ultimately lead to the classical equations of quasigeostrophy. In particular if 
C ^ 1 then this limiting form identically leads to classical QG. We see that the first non-trivial inclusion of magnetic 
effects comes in when C is an order 1 quantity. In this event, the resulting equations will be called magnetoquasigeostrophy 
(MGQ). The completed derivation of the latter is found in Section 01 



3.2.2. Magnetostrophic balance 

When magnetic effects are relatively strong there exists another more general lowest order balance between Coriolis, 
horizontal pressure gradients and the horizontal components of the Lorentz force. This magnetostrophic balance has been 
explored in the context of the Earth's geodynamo (Acheson & Hide 1973, Fearn 1997, Moffatt 2008, to name just a 
few). Inspection of ([TO)) shows that in order to bring the Lorentz term in on the same order as the Coriolis term requires 
C = 0(Ro-^). Thus, we define a new number 

A 



which we shall, henceforth, call the Acheson Number. For this three-way balance to occur the Acheson Number must 
be an order 1 quantity. Similar inspection of (jl9|) shows that in order for the horizontal pressure gradient to enter the 
mix we must have the combination (5-Bu^ — ©(Ro). Inspection also shows that if the resulting scaling arguments are 
consistent with the resulting equations 6 must be O (1) (see below). Thus we define a modified Burger number Bu such 
that Bu = Ro^/^Bu, where 



Bu 



With the assumption that (i) all quantities may be expanded in powers of Ro (as we did in the previous section) and (ii) 
the lowest order horizontal velocity and magnetic field components are z-independent in a layer, then we find that the 
magnetostrophic balances in the horizontal directions are 

2z X uf^ = -BVlvnf ) + -(bf ) . V)hf\ (21) 
Pi Pi 

The horizontal divergence of the lowest order velocity is not zero in general. In fact. 



^.4" + d.vf^ = ^ {d.b^ + d.b^ + b|") . V) ^ 0, (22) 



in which J^^"^ is the lowest order vertical current in the layer, J^^^ = dxb^y^ — dyb^^ . Unlike what we encounter in the MGQ 
scalings, it follows that the lowest order vertical velocity cannot be zero at lowest order. Since S measures the vertical 
variations of each layer height and since the latter is related to the vertical velocity, in order for the magnetostrophic 
scalings to be consistent and to have the lowest order horizontal divergences to be non-zero it must be that S = O (1), 
which is to say vertical motions are driven by vertical currents due to the Lorentz force at lowest order. The continued 
analysis of the magnetostrophic equations and development of a self-contained equation set is presented in Section [S] 

In Table [T] we summarize the relative scalings of the relevant quantities that produce the reduced equations discussed 
to now (i.e., QG, MQG and magnetostrophy). 

3.3. Analysis of the incompressibility and induction equations 

The scalings assumed at the beginning of this section result in the incompressiblity equation in each layer rewritten as 
dxUi + dyVi + dzWi = 0. (23) 
The induction equations are also similarly written as 

{dt + u,-V + w,d,)h, = (b, • V + b,,d,)n, (24) 
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with the vertical field component bzi relating to the other field components through the divergence free condition 
dxbxi + dybyi + dzb^i = 0. (25) 

4. Detailed derivation of the equations of magnetoquasigeostrophy 

We formally develop the scaling analysis and equation reduction that we began in Section [3.2.1l Quantities in each layer 
i are expanded in the following way, 

u, = uf)+RouW + ... 
b, = bf +RobW + ... 

n, = ' + Ro np' + • • • , (26) 

where the superscript denotes which order of the Ro expansion the term represents. As we stated before, insertion of 
these expansions into (|19p results in ([20)1 to lowest order. Since we are considering a two-layer problem we have by layer 
that 

= -a,|^nl°), «(°).a.|^nr, (27) 

Zpi Zpi 



and 



.f = -a,|^iif .f^a.f^nf, (28) 



It will prove to be more convenient to write the pressure fields 11^ in terms of a streamfunction ipi such that 

(29) 



(0) ^ Bu^^(o) ^ Bu^ r 
' 2pi 1 2 



Pi " V Pi 



and 



-l^nf .^/.(°), (30) 

where the height deviations have been similarly expanded in powers of Ro, i.e. 
/i, = /jf)+Ro/.(°) + ... 

This rewriting expresses the relationships more transparently 

.r = -a,4°\ = = vi'' = dA''- (31) 

Evidently the fields u|°' are 2-dimensional incompressible flow independent of the vertical coordinate z. This means from 
the incompressibility equation (j23p 

we get 

d,wf^ = 0. (32) 
Thus the leading order vertical velocity term is zero motivating its expansion to be 

w, = Ro wf^ + • • • (33) 

Note that this is necessarily consistent with the assumption made at the outset that 5 — O (Ro) . In other words this says 
that if the vertical velocities are small then so must be the vertical variations of height. 

We turn to the leading order equation for the horizontal compoenents of the magnetic field and we find 

(a. + uf'.v)bf = (bf .v + 5i°)9.)uf). 

However, since the horizontal velocity components are independent of z the above equation reduces to 
(9. + uf).v)bf).bf .Vuf), (34) 
if the totality of all the layers are not moving. 
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which implies that hf'^ is independent of z as weh. The leading order expansion of the source free condition ([25] 



d.l^+dyl^^~d.l^ (35) 

relates the leading order vertical field 6^'^'' to hf*^ . We note that it means 6^°' is at most a linear function of the coordinate 
z within each layer. Furthermore, and quite unlike the geostrophic balance situation we encountered before for the leading 

order horizontal velocities, the horizontal magnetic fields are not strictly two-dimensional, i.e. b^z'^ is not identically zero 
(Oilman, 2000). It follows from this that there are modes of dynamical activity which will allow the divergence of the 
leading order horizontal magnetic field resulting in a certain amount of "breathing" into the vertical field component. 
At order Ro of the incompressibility condition we find 

+ = . (36) 

The remaining order Ro terms of the horizontal momentum balance equation pop are 

(9, + uf ■ V)uf ) + 2z X uW + /3yz X uf ) = -Bu^lvft^ + C(bf ' • V)bf (37) 

Pi 

where we have explicitly used the fact that b^°^ is independent of z in writing (j37p . Because an examination of the next 

order expansion of the vertical momentum balance equation reveals that n|^'' is independent of the vertical coordinate, 

it follows that the next order horizontal velocity corrections u^^-* are also independent of z. We may therefore take the 
curl of the above equation by operating the y-momentum component by dx and subtracting from it the dy operation 
upon the x-momentum equation to get 

(a. + uf).v)(a..f -a,.r)+2f^«. 

C,(bf ) . V)(5.6j,°) - dyb^^) + a{dj^^ + dyhf){dxb'^S ~ 9yb^h- (38) 
with Ci = C I Pi and where, for notational convenience, we have defined 

= -{djP +dyV^^) (39) 

In order to proceed we must determine how f^f ^ relate to hf\ To do this we begin by determining the equation for 
each of the level heights hi. The previous scaling analysis leads to the non-dimensionalized version for them which, at 
leading order, reveals that 

{d^ + nfyV)hf ^w'^\z = H,) (40) 

where care must be taken since one may evaluate the motion of the j^^ surface using the velocity fields from layer i. 
Specifically we explicitly write out the solution to the vertical velocities for each layer. In the bottom layer we have 

zi-^'^ = (41) 
in which we have implicitly forced the vertical velocity to be zero at the bottom ^ = 0. To leading order it follows that 

wf\z ^ Hi) = Hwn^^\ (42) 
When viewed from the lower layer then we have the motion of the lower layer's interface to be given by 

(d, + nf^-v)hf^^H,,n^^\ (43) 

For the upper layer we have 

= u;W + f^Wz, (44) 

where Wjo'* * constant velocity to be determined below. Now it is evident that the motion of the bottom boundary 
when viewed from the upper layer is 

(dt + u(°' • V) hf^ = ^4;,) + i/ior!^'^ , (45) 
while the motion of the upper boundary when viewed from the upper layer is 

(dt + • V) = + ^20^!^'^ . (46) 
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Subtracting these two expresses the evohition of the width of the upper layer in terms of ri2^\ that is to say, 

(d, + uf^ ■ V) (4°^ ^ ) = (i?20 - H,o)4'\ (47) 

which, we note, is purely in terms of the velocity fields U2'^''. We are now in a position to rewrite (j38|) layer by layer by 

replacing all instances of 
begin with the bottom layer 



replacing all instances of with its appropriate expressions in terms of given by (|43)) and ((47|) . To this end, we 



while for the top layer we have 



2h 



(0) 



H 



10 



(48) 



(dt + u! 



(0) 



• V 



(0) r,^2 ~ "l 



(49) 



In order to simplify the expressions appearing above we rewrite the heights hf''^ in terms of the streamfunctions 

2 



,(0) 



2 /V^r^-tV^f 



,(0) 



Bu^ 



1 _ PI 
pi 



Bu 



(0) 



We may now rewrite the equations in more transparent form 

(a. + ur • v) gr + - c (b^ ■ V + ^\) 4'\ 



c 



P2 



where the potential vorticity in each layer is defined by 



pi 



^10 



l(0) - T72 



in which the Lio and L21 denote the Rossby radius of deformation for layers 1 and 2 respectively, 



1 



1 



^?o Bu^ifio (1 - ) ' ^21 Bu2(F20 - Hio) (1 - 



(50) 

(51) 
(52) 

(53) 
(54) 

(55) 



The form presented above diverges slightly from the form presented in Vallis wherein the upper layer is constrained by 
an upper boundary while here we allow for the upper layer to undulate freely. The difference then is that the ratio P2I P\ 

appearing on the RHS of defining Q^^^ would be replaced by 1. We have also defined the currents jf^ appropriate 
to layer i 



•J\ — i^xUyi iJyUxl ' "^2 — ^xUy2 ^y"x2 ' 

and the corresponding layer densities g^j] , 

^(0) ^ Q dO) Q AO) JO) ^ , (0) Q AO) 



(56) 



(57) 



All the terms in these equations, as appearing in the summary Section 6, will be expressed with their individual super- 
scripts "(0)" removed. Additionally, the expression for Bu is rewritten in terms the definition of the Rossby Deformation 
Radii. 
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5. Further development of the magnetostrophic limit and the derivation of a closed set of equations 

The leading order magnetostrophic balance (PT|) leads to a diagnostic specification of the horizontal velocities which we 
explitly write here 

2.^' = - (^dM' + -bf^ ■ Vl^ 1 , (58) 




2ur = :i^a,ri|n^bf).vC. (59) 
Pi Pi 



We shall assume henceforth that the horizontal components of the velocity and magnetic fields are z-independent within 
a given layer. Thus, the leading order induction equation is the same as ([M]) which we rewrite here for convenience, 



The horizontal velocity components are expressed in terms of b^*^^ and the perturbation height fluctuations hi are 

contained in the perturbation pressure fields tl\^\ Because the horizontal divergence of the lowest order velocity field 
was shown to be non-zero and subsequently assumed to be independent of the vertical coordinate z, we can write the 
leading order vertical velocity to be 

where, as before, w^-g-* is a z-independent undetermined vertical velocity. Because of the relationships set up between the 
velocities according to the leading order continuity equation 



together with the help of (j22p . we can express the quantity olf^ in terms of the horizontal magnetic field quantities as 

^-^^(dJ^+dJS^X^-v)jr. (60) 

We develop an equation for the height fluctuations in the same manner as was done in the development of the MGQ 
equations. Because we have argued the fluctuations are on the same scale as the mean height of any given layer (i.e. 
because 5 — 0(1)), instead of expressing the height Hi as the sum Hio + Shi we shall simply stick with the expression Hi 
since it now makes no difference as both quantities are the same order of magnitude. Thus the vertical motion of layer 
height with position z = Hi as viewed from the layer with index i is to leading order 

(a. + uf).v)i/. = ^|^^^^=^(°) + i/Ar. 

Similarly the motion of the layer height with position z = as viewed from the same layer with index i is to leading 

order 

(a. + uf ' . v) ^ = -S^ + ^-i^^r^ 

Subtracting these two equations and reordering the results reveals 

{d, + uf) . V) In (H^ -H„^,) = = - A (a^foCo) + a^,(0) + ) • v) J^. (61) 

As presented, this magnetostrophic limit involves evolving equations (|22p and (|6ip with the horizontal flow quantities 

determined by the magnetstrophic balances explicitly given in (1581) and (|59p . where the perturbation pressure fields n^'^'', 
given in ^TE\\. are re-expressed in terms of Hi. We remind the reader that (IT51) explicitly shows the pressure field quantities 
for a two- layer system. In general, for layers i < N with hi expressed in terms of Hi, we would have, 

N 

nf'>^p,H,+ J2 Pk{Hk-Hk-i), (62) 

k=i+l 

and for layer N we have n^"* = pnH^. 

To illustrate how these equations appear we write them out for a single layer atmosphere where we drop the super- 
scripts, 

((9t-l-ui • V)bi =bi • Vui, (63) 
{dt + ui . V) Inffi = -(A/2) {d^b^i + d^byi bi • V) Ji, (64) 
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in which the velocities and currents are diagnosticahy given by 

2wi = -Bud^Hi - Abi • \/b^, (65) 

2ui = BudyHi + Abi • V6.y, (66) 
Ji = dxby- dyhx. (67) 

Note that these equations are governed by the two parameters, the modified Burger number Bu and the Acheson number 
A. 

6. Summary and brief discussion of the equations of magnetoquasigeostrophy 

The result of the lengthly procedure detailed in Section 2] are summarized for two layers. The equations for the potential 
vorticity in each layer are given explicitly, 

{dt + ui • V) Qi + Pv^ = ifiiC (gj"^ + bi • v) Ji, (68) 

{dt + U2 • V) Q2 + PV2 = f2C (^^^ (g^™^ + b2 • V) J2, (69) 

in which the horizontal velocities in vector form are given by = Uik + viy. The gradient operator V is now understood 
to be two-dimensional so that, for instance, • V = Uidx + Vidy. The potential vorticities for each layer, Q^, are given by 

72 J. 1 (^,. P2_ 



V^^i - — - , (70) 

^10 V Pi 



Q2 - VV2-7^(V^2-V'i), (71) 



1 

^21 

where the streamfunctions ■0^ relate to the velocity in each layer according to 

Ui = -dyipi, Vi = dxipi, (72) 

and where the 2 dimensional Laplacian operator is given by — > d^^^y. Note that because the leading order horizontal 
velocities are geostrophic it follows that their horizontal divergences are zero, i.e. 

V • u, 0. 

The (non-dimensionalized) Rossby radii of deformation for each layer are given by 

.2 _ .g^glO P2\ ,2 _ 9'H{H2G -Hio) P2\ 

where pi and p2 are the densities of each corresponding layer. The dimensional value of each level height T-li are given 
respectively according to their non-dimensional fractional measures HiQ. In other words if Tiio represents the vertical 
coordinate where the transition from density pi to pi+i occurs in steady state then the following may be defined Hiq = 
TiiofH, H20 = 'H2o/'H. Since only two layers are considered here the value of H20 = 1 as H20 = 'H by construction. It 
follows that Hio is some number less than 1. The horizontal magnetic field in each layer is denoted in vector form with 
bi = bxii- + byiY and the equations for their evolution are given by 

(at + ui • V)bi = (bi • V)ui, (74) 
(at + U2-V)b2 - (b2-V)u2. (75) 

As in the study considered by Oilman (2000), and unlike what is encountered for the horizontal velocity components, the 
horizontal magnetic field components are not automatically divergence-free. Of course, the source free condition is met in 
three dimensions so that a non-zero divergence of the horizontal magnetic field components will result in the generation 
of a linearly varying (with respect to the vertical coordinate) vertical magnetic field bzi, as expressed in relationship 
l(55|) . It so happens that for these orderings the vertical field (if there is one dynamically, see below) effects the dynamics 

diagnosticahy through the pseudo-source term q^"^^ defined by 

gl") = V • bi = dxbxi + dybyi, (76) 

and its influence can be seen in the Lorentz terms in the equations for the potential vorticity (|68II69I) . The source term 
q{m) -j-j^g pseudo-magnetic monopole distribution referred to in the Introduction. Finally the z-directed currents for 
each layer are Ji and defined by 

Ji = dxbyi - dybxi- (77) 
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An equation for the source terms qf^^ may be developed by operating each equation (j74p and (|75p by the divergence 
operator and making use of the fact V • — 0. This procedure resuhs in, 

(9t + u,.V)<z(™) =0. (78) 

This equation says something important: that if qf^^ {x, y) — initially then it remains zero for subsequent times. When 
the case, the consequence of this is that the horizontal magnetic field components obey a divergence-free condition and 

the evolution of the magnetic field is strictly two-dimensional derivable from a fiux function. Note also that qf^\x, y) — 
solutions of the MGQ equation set form a self-contained subspace of all solutions of the MGQ equations. When 
examination is focused wholly on this subclass, i.e. when 

gM=g(™Uo, (79) 

it automatically follows that V • = 0. Consequently, the flux functions relating to the magnetic field components are., 

hxi = -~dy<pi, byi = dx(t>i, = -9j^02, by2 = dx(l)2- (80) 

Instead of solving the vector equations (fM)) and ([75)) . the following simpler scalar equations must solved, 

(dt + ui • V)0i = 0, (81) 
{dt + U2 • V)02 = 0. (82) 

In terms of fiux functions, then, the currents in each layer are given by 

Ji - V^^i, J2 = V202. (83) 

Thus the evolution equations given by (|68l69p . with gl™-* set to zero, and (|81I82I) . together with all of the supporting 
ancillary definitions form the basis of the equations appropriate for this subclass. These are the same equations analyzed 
by Oilman in 1967-8. 

In an upcoming series of publications, we study the linear and nonlinear response of these equations which are 
amenable to relatively unencumbered theoretical analyses using well-known techniques used in both the meterological 
and fluid dynamics literature for studies of the QG equations. To this end, we shall use as a guide and expand upon the 
work laid out by Oilman and collaborators. Our initial investigations will focus upon analyzing the physical stability of the 
compact magnetic vortices observed to emerge in the simulations of Oho (2008) and to hopefully, as a result, understand 
the nature of those structures and answer why they may be manifesting in those simulations with the observed robustness 
reported. 
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